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A variety of open quantum networks are currently under intense examination to model energy 
transport in photosynthetic systems. Here we study the coherent transfer of a quantum excitation 
over a network incoherently coupled with a structured and small environment that effectively models 
the photosynthetic reaction center. Our goal is to distill a few basic, possibly universal, mechanisms 
or "effects" that are featured in simple energy-transfer models. In particular, we identify three 
different phenomena: the congestion effect, the asymptotic unitarity and the staircase effects. We 
begin with few-site models, in which these effects can be fully understood, and then proceed to 
study more complex networks similar to those employed to model energy transfer in light-harvesting 
complexes. Our numerical studies on such networks seem to suggest that some of the effects observed 
in simple networks may be of relevance for biological systems, or artificial analogues of them as well. 

PACS numbers: 



I. INTRODUCTION 

The transport of electronic excitations over biological 
networks of chromophores is the relevant mechanism for 
the light-harvesting step of photosynthesiJ^^^. Recently, 
long-lived quantum coherent oscillations have been ob- 
served in ultrafast experiments carried out on several bi- 
ological systems, even at room temperatur^^EH^ One of 
the key features of these exciton-transfer networks is their 
open nature, namely, that their coupling with the protein 
vibrational environment is, arguably, the dominant effec- 
tor of transport in these systems. The interplay of uni- 
tary dynamics and the system-bath interaction has been 
predicted to be beneficial to the network functionality at 
biological conditionJ^^HH] xhe different mechanisms that 
lead to this environment- assisted quantum transporf^ 
are still under vigorous exploratioiP^. Realistic numeri- 
cal modeling of these open quantum networks is, to some 
extent, possible and currentl y act ively pursued in the 
physical chemistry communitjJ^nillHlil, Nevertheless, the 
physical chemistry and quantum information co mmunity 
has learned much from simple Markovian models^SIMEoi 

In this paper, motivated by the above, we will in- 
vestigate a few simple yet illuminating models of open 
quantum networks in order to identify a handful of basic 
mechanisms or effects that are featured in fully analyz- 
able toy models and that may persist for larger, more 
complex quantum transport networks. In particular, we 
will focus on coherently-coupled qubits subject to dissi- 
pation/dephasing and irreversibly connected to an auxil- 
iary quantum system. The role of this latter is to model 
the reaction center of light-harvesting complexes, where 
the electronic excitation is separated into an electron and 
a hole and the charge-transfer stage of photosynthesis 
begins. Of interest to us is the reaction cen ter o f the 
LHl-RC complexes present in purple bacterieP^ll. We 
will adopt a Markovian master equation of the Lindblad 
form to describe the overall system dynamics. Different 
energies, or equivalently time-scales, will enter the def- 



inition of the Liouville superoperator C. The interplay 
of these time-scales controls the non-trivial phenomenol- 
ogy that we explore in this manuscript. Finally, singling 
out a few intriguing, possibly universal features of such 
a phenomenological landscape is the goal of the simple 
calculations presented in this paper. 

In the next three sections (|llj III and IV) we will 
consider different toy models consisting of few sites or 
chromophores (modeled as quantum two-level systems, 
or qubits), manifesting particular features which can be 
fully understood by analytical calculations. See Fig.[l]for 
a cartoon picture of the various networks considered. In 
section [Vj we will consider more realistic networks bor- 
rowed from models of light-harvesting complexes. Via 
numerical simulations we will show that these effects may 
persist in more realistic systems. 



II. THE CONGESTION EFFECT 

In exciton and electron transfer events, there can be de- 
lays in energy transport due to the timescales of the bio- 
logical process. A particular element might be shut down 
while transport takes place, effectively making an exci- 
ton or electron wait until the transport is possibl^^. In 
the following section, we will describe this phenomenon 
in model systems and characterize it as the congestion 
effect. 

In the standard modeling of incoherent (and irre- 
versible) transfer of excitations from one site to an- 
other, the Forster electromagnetic coupling mechanism 
permits the transfer of populations at a given rate 7. 
If the dynamics is described using a Lindblad form 
p = Cl (p) , where Cl (X) = LXL^ - {L^L,X] /2, this 
can be accounted for by a jump operator of the form 
L = s/^a^ ® (7+, where a* are Pauli ladder operators. 
For our purposes, such a Lindblad description is phe- 
nomenological. Site 2 could model, for example the re- 
action center of LH-II described above. In this section. 
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L = ^ya^ o-j L=4yo-^ 

Figure 1: Summary of the toy-networks analyzed analytically 
in sections Ull and Hill 



we explore possible congestion effects tfiat arise from the 
dependence of the transfer rate on the number of exci- 
tations involved, in the same way traffic flow might be 
inversely proportional to the number of vehicles present 
on roads. 

Incoherent transfer /; • ». Before turning to ana- 
lyze the possible implementations and consequences of 
such an effect, let us summarize the Lindblad opera- 
tors for incoherent Forster transfer among two sites, 
L = G2 ■ This process can be pictorially visualized 

by the following diagram, • • (see also Fig. The 
quantum master equation is given simply hy p = Cl (p). 
We denote by n the population operator satisfying nj?/) = 
r]\r]) with 77 = 0, 1, and by n its possibly time-dependent 
expectation value for excitations, i.e. n = (n) := tr(np). 
Since the effect of the Lindbladian is to transfer a par- 
ticle from site 1 to site 2, the total number operator 
is a conserved quantity. We therefore obtain a differ- 
ential equation for the population in the following way: 
first note that rii = tr (riip) = tr [ni^L (p)] ■ Given that 
ni + n2 = ntot is constant in time, it suffices to analyze 
the population of site 1, ni = —jui -f 7(nin2). Now note 
that in the single-particle sector, titot = 1, (^1^2) = 
(to see this use n^^^ = ntot + 2nin2), leading to a trans- 
port equation hi — —7711 that can be readily solved 
for the population at sites 1, ni (t) — e~^*ni (0) and 
2, 712 (t) = n2 (0) + (1 — e~^*) ni (0). The jump operator 
achieves precisely what we expected: the population in 
site one decreases exponentially at a rate 7 and the pop- 
ulation of site 2 increases accordingly. The same result 
could have been obtained by solving the (16 dimensional) 
differential equation for the full density matrix. Starting 
at time zero with p (0) — {pij} the time-evolved density 
matrix p (t) is 
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It is interesting to note that for some entangled initial 
states the asymptotic density matrix p (f — )■ 00) is still 
entangled. The process Cl cannot, however, create en- 
tanglement. 

Incoherent transfer II: • □. To model the con- 
gestion in the reaction center, let us now substitute the 
second qubit with a larger 2s + 1 dimensional space. 

For this case, we can model a particle conserving trans- 
fer process with a jump operator given by L = y^cr^ S2 
where S2 is a raising operator of the irreducible spin 
s representation of SU{2). The population at site 2 
is D^2 = 82+ si. Once again, since the total particle 
number ntot = ni + 9I2 is conserved in a given particle 
sector, one has ntot (t) — "-tot- We then obtain the fol- 
lowing differential equation for population at site 1: hi = 
-7(niS2"S'2^). By noting that 5-^52^ = (fn2 + l)(2s-9l2), 
and employing OT2 = "tot-ni, = 7i2^t-2"-tot-l+2ni, 
and n^ = ni, we obtain an explicit differential equation 
for ni: 

ni = -7'^tot [(2s + 1) - ntot] ni 
ni + N2 = ntot- 

Excitation transfer now occurs at an effective rate 
which depends on the total population; 7off = 
77T,tot [(2s + 1) — ntot]. Note that < ntot < 2s + 1 and, 
correctly, 7cff (ntot = 0) = 7cff (ntot = 2s + 1) = 0, i.e. no 
transfer takes place when the network is either com- 
pletely empty or completely full. The maximum transfer 
rate is attained when the condition ntot — (2s + 1) /2 
is satisfied. The lesson we get from this slightly modi- 
fied example, is that transferring excitations to an object 
with more than just two levels, is likely to result in a 
population dependent transfer rate. 

Interplay between coherent hopping and transfer: • -H- 
• • » . We will further illustrate the concept above 
by considering a variation on the theme. We consider a 
coherent-hopping Hamiltonian on four sites of the form 
H — (J/2) (c^^cr^ + h.c.) that acts on the first two sites. 
The excitations are transferred irreversibly from site 2 to 
site 3 via a quantum jump operator L = -y/7cr^(T;^ and 

subsequently from site 3 to site 4 with = ^/ib'^-i^^t- 
J is the coherent coupling strength. In the following, we 
explore the interplay between the two incoherent transfer 
rates 7 and 7f,. Let us focus on the population at site 3, 
n^{t). The effect of 7^, is that of removing excitation 
population from site 3. However when becomes large, 
excitations are rapidly transferred to site 4 inhibiting the 
effect of Lfc {Cb (p) — !• 0). This results in a non-trivial 
non-monotonic effect as a function of jb- This feature 
can be visible only if we have at least two particles in 
the network. Let us then consider the following initial 
(pure) state with excitations localized at sites 1 and 2: 
jl, 1,0,0). As shown in Figure 3, in this case, the time- 
evolution of the populations takes the following form: 
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Figure 2: Top panel: The population 713 for the case described 
in section [n], as a function of time and jt- In this model, the 
initial state has two excitations at sites 1 and 2: j 1,1, 0,0). 
The parameters for the model are J — 1, — 0.1. Bottom: 
panel: Slices of the same plot at different times are shown. 
The non-monotonic behavior of the population as a function 
of the rate 76 is evident at small values of it. 
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where Ci, C[, Bi are only functions of J, 7, 7;,, and 
C3, C3, Bi are functions of time as well. Finally a; = 
^ ^2 _ 4 j2 ^ resulting in an imaginary eigenvalue of the 
Liouvillian for 2\J\ > 7. This in turn shows up in an os- 
cillating behavior of the populations as a function of time. 
In Figure [2] the behavior of population 3 as a function of 
time and 7b is plotted for the given values of J and 7. For 
large values of t, one can observe a non-monotonic behav- 
ior as a function of 7f, emphasized in the bottom panel of 
Figure [2] This behavior can be qualitatively understood 
as follows. Consider the behavior of as a function of 
7b for a large fixed time t. Since the effect of 7b is that of 
taking away particles from site 3, first decreases when 
7b is increased from zero at fixed t. Anyway, if 7b is fur- 
ther increased, excitations are taken away at a faster rate 
and transferred to site 4. This means that at the fixed 
time i site 4 tends to get full for large 7b, thus inhibiting 
the effect of Lb- Population 713 then increases with 7b. 



Figure 3: ns, as a function of time and 76 . Initial state has one 
excitation at sites 1: jl, 0,0,0). Parameters are J = 1, 7 = 
0.05. 



When 7b is further increased, site 4 becomes effectively 
full and Lb is turned off, the population becomes then 
independent of 7b and 0.3 saturates. 

For the sake of completeness we also consider the solu- 
tion with one excitation localized at site 1, i.e. |1, 0,0,0) 
at time t ~ 0. In this case the time-dependence of the 
populations is. 
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where Ai are time independent functions of the param- 
eters. One can see in Fig. |3] that the non-monotonic 
behavior of 713 as a function of 7b is for this initial con- 
dition absent. As expected, since in the network there 
are no-excitations enough to fill the reaction centre, the 
"congestion effect" is now absent. 



III. THE STAIRCASE EFFECT 

In this section, we explore the situation where excitons 
are fed into a quantum network at a given constant rate 
7in and are extracted at a rate 7out- 

This model can be justified by the fact that some 
photosynthetic complexes such as purple bacteria and 
green-sulfur bacterigP^l live in low-light conditions. The 
electron-transfer event that occurs in the reaction center 
is a process that takes place in the order of picoseconds. 
We therefore take the common practice of modeling the 
reaction center as an incoherent trap-'^. 
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explicit expression of the asymptotic state is 



Figure 4: J = 2, 7in — 0.2 (injection), and 7out = 0.3 (dissipa- 
tion). Bottom: parametric plot with same parameters. The 
red dots correspond to times given by T„ = (l/2-|-n)ro, 
To — 2n/ijj, (n = 0, 1, . . .) and the correct frequency is 

u) = \J 4J-^ — (7in — 7out)^. As explained in the text using 
general arguments, To = O 



Injection-extraction: • o • Here, we con- 

sider the simplest model for the injection and extrac- 
tion of an exciton. The model corresponds to two sites 
coupled coherently via the hopping Hamiltonian, H ~ 
{J/2) {^a^ a2 + crt(J2) ■ Besides the coherent evolution 
term, an incoherent injection of excitons is given by a 
jump operator Lin = y/jinO'^ which injects particles at a 
rate 7i,i and a corresponding incoherent extraction term 

Lout — VToutCTj . 



The corresponding 16 x 16 Lindblad superoperator ma- 
trix can be diagonalized. A complex eigenvalue with a 
non-zero imaginary part gives rise to oscillating behavior 
in the populations when |7iii — 7out| < 2 | J|. 

Let us first concentrate on the asymptotic state of the 
evolution, p(t ^ oo). Solving £tot if) = 0, one realizes 
that the asymptotic state is unique and independent of 
the initial state. Although this feature is expected in nat- 
ural physical systems and follows, for instance, from the 
detailed balance hypothesis, it is not necessarily satisfied 
in our simple toy models (see e.g. Sec. [TT|). 

In the standard basis, {|1,1), |1,0), |0, 1), |0,0)}, the 
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The only non- vanishing correlations are (cf (t|), (trf ) and 
(erf). Thus this state is separable but has non vanishing 
classical correlations: (cr^crl) — (o'J)(cr|) ^ 0. Equiva- 
lently, the asymptotic state is a classical mixture of states 
with definite populations. 

Having p (oo) we can compute the asymptotic popula- 
tions: 

f ^ 7in {J^ + 7in7out + 7out) x 
ni[oo) = — — , , ^„ r (1) 



n2 (oo) 



(7in + 7out) [J^ + 7in7out) 

7inJ^ 

(7in + 7out) [J^ + 7in7out) ' 



(2) 



A few simple facts can be directly seen from equa- 
tions ([T]), ([2|. First, for small 7in populations deviate 
by O (Tin) from zero; vice versa for 7out small popu- 
lations deviate by O(7out) from one. Instead, when 
J is small excitations get loaded at site 1 but take a 
long time to reach site 2 so that ni = 1 — 0{J^), 
n2 ~ O (J^). Finally, for very large J both populations 

tend to ni ~ ^2 = 7in/ (7in + 7out) + O (J~^). 

Let us now turn to the dynamics and consider first the 
most interesting case namely when the initial state is the 
empty state |0,0). A typical result is shown in figure |4] 
An interesting feature clearly emerges: when population 
ni increases, n2 stays almost constant and vice-versa. 
Such a feature is particularly evident in the paramet- 
ric plot. In the lower panel of figure [4] we also stressed 
another peculiarity of this process: the time needed to 
increase a given population when the other is constant 
(i.e. the horizontal and vertical steps between two red 
dots in Fig. |4]), is always the same. We call Tg this new, 
emerging, time-scale. The description of the entire pro- 
cess then is the following. First particles are injected at 
site 1 and population at site 2 stays zero until a time 
To/2. Next, for To/2 < t < 3/2To the situation is re- 
verted and population 2 increases while population 1 re- 
mains constant. The process continues in this fashion 
until an asymptotic state is reached. Given the shape of 
the curve in Fig.[4]we refer to this situation as "staircase 
effect" . The emerging time-scale can be given a physical 
interpretation considering the limit when both injection 
and extraction rates are very small. In this case the only 
time-scale of the system is given by the time needed for 
the excitations to hop from site 1 to site 2. This time 
is given by Tq « AE~^ = 0(J~^). In general, if we 
substitute the two sites with an open chain of length L, 
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asymptotic populations are given by 
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Figure 5: |1,0). Parameters are J = 2 7in — 0.2 (pump) 
and 7out ~ 0.3 (dissipation). Below; parametric plot, same 
parameters. 



using the same argument we expect (at least for small 
7in, Tout) that Tq will be the time needed for the exci- 
tations to travel from one side of the chain to the other, 
i.e. To ~ L/v where v is the velocity of quasiparticles. 
Of course this picture can be correct only as long as a 
quasi-particle description applies (cf. Sec. VCl. 



Let us now consider the injection-extraction dynamics 
with an initial state |1,0), i.e. at time zero the injection 
site is occupied. A typical (in the oscillating regime) sce- 
nario is shown in Fig. [5j Starting with an initial state 
|0, 1), the situation is almost identical with ni and n2 
interchanged. In fact, one can show that for initial states 
with one definite excitation, populations at any time sat- 
isfy the following duality relation 



ni (7in, 7out) = 1-^2 (7out, 7in) • 



As previously explained the asymptotic populations do 
not depend on the initial populations and are still given 
by equations ([I]) and The parametric plot in the 

lower panel of Fig. [5] shows that with this initial condition 
the staircase effect is not present. 

Three-site injection-extraction: ^ • o • <4 • A 
slight generalization of the above idea is given by a three 
site chain with injection on the first site and extraction 
on the third. For simplicity we consider a uniform chain 
with equal couplings J12 = J23 = J- In this case the 
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Note that populations ni and are the same as ni, 
712 in the previous, two-site case. Starting from the to- 
tally empty state, the asymptotic state is reached in a 
similar manner as in the two-site case. In particular, 
the parametric plot of the injection and extraction sites 
(ni (t) , {t)) displays a staircase shape exactly as in the 
two site case. As we will show in Sec. |V C[ this feature 
survives even in a longer chain, and is to some extent re- 
sistant to small static diagonal disorder and dephasing. 



IV. ASYMPTOTIC UNITARITY 

Another effect we want to study is the possibility that 
a coherent dynamics (or sub-dynamics) may emerge out 
of a dissipative or partly incoherent dynamics. To make 
things more clear let us immediately discuss the simplest 
example showing this feature. 

J 7 

Hopping and transfer: • o • •. The model 
consists of three sites (qubits). On the first two 
sites acts a coherent hopping of the form H — 



(J/2) (crfcr^ -f crj*'cr^). On top of that, particles are 
transferred irreversibly from site 2 to site 3 via a jump 
operator given by L = \/l'^2'^t- clear that, if a 

particle sits at site 3 the incoherent part of the dynam- 
ics is not effective, that is £l [pi2 — 0- If we 
start with an initial state |1, 1,0) with sites 1 and 2 oc- 
cupied and site 3 empty, for effect of the dynamics, site 
3 will get populated at a rate 7, and on the first two 
sites there will remain one particle coherently hopping 
back and forth. By this we mean that for a sufficiently 
large time the evolved state will be similar to a coher- 
ent evolution: p{t) = e*^'°' [p] ~ e~**-^pe'*^ =: p(t). 
For what concerns the state p we only know that it will 
contain one particle; it can be obtained by evolving back 
unitarily p (t), i.e. 



P ■ 



lim 



itH 



p{t), 



itH 



Indeed, if the dynamics becomes unitary, the above limit 
is well defined. Notice that p is nothing but the station- 
ary solution of the original master equation in the interac- 
tion picture associated with H. The same reasoning can 
be done for the subsystem consisting on sites 1 and 2, 
i.e. we can define pi^2 by evolving back unitarily pi^2 (t)- 
Since H does not act on site 3 we have pi^2 = trsp. An 
explicit computation confirms that p = pi 2'8'|1)(1|, i.e. in 
the equivalent, unitary dynamics, one particle sits at site 
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3. The explicit form of pi.2 in the standard basis is 
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This state is a quantum superposition of one-particle 
states with m = 1/2 + 7^/2 (j^ + 7^) and ^2 = 1/2 - 

What are the possible indicators of asymptotic unitar- 
ity? Since the purity is constant under unitary evolution, 
one possibility is to look at the purity of the total sys- 
tem or of some part of it. The time-derivative of such a 
quantity will then be close to zero, for approximate uni- 
tary evolution. Since for Lindbladian evolution the purity 
derivative is 9ttr (p^) = 2tr [p (t) £tot (p)], this definition 
has the advantage of being numerically stable. In our toy 
model we have 




Figure 6: {ip \pi,2 (t) for the model considered in the 
text. Parameters are J — 2, j — 1. 
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2(J2 



-+ 



-26-^* ( J2 +72+72 COS (Jt)) + e-27t (3 J2 + 4^2) 
2(J2+72) ■ 

Unfortunately the purity tends to a constant whenever 
the solution tends to a constant, as happens, for instance, 
along the natural process reaching the asymptotic state. 
In other words, the smallness of the purity derivative is 
a sufficient but not necessary condition for asymptotic 
unitarity. 

Another possibility is to measure some distance be- 
tween the actual state and the one obtained with unitary 
evolution: \\p{t) — p{t)\\. Once again, we might as well 
restrict to a particular subsystem. Using the operator 
norm the result for our toy-model is particularly simple 
and illuminating 



This confirms our initial intuition: the dynamics becomes 
unitary at a rate 7. This approach has a very clear mean- 
ing but has the disadvantage of being computationally 
demanding as it requires the computation of a matrix 
norm and the evaluation of p{t). A simpler alternative 
is the following. 

Consider the spectral representation of the Hamilto- 
nian H = ^n^n\n){n\. If the total evolution becomes 
similar to a unitary evolution, the matrix elements of the 
density matrix in the eigenbasis \n) evolve in time like 
phases: 



{n\p (t) \m) ~ (n|p (t) \m) = e 



-it{E„-Err,} 



n\p\m) 



In our model the eigenbasis of the two-site Hamiltonian is 
{|0,0), |1,1), IV'^^) = (|1,0) ± |0,1)) /V2}. For instance, 
one can show that 
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Pictorially the parametric plot of the real and imaginary 
part of this matrix element folds on a circle (of radius 
l/y/J^+l^) after a time 7 ^ (see Fig. [6|. 

This method to mark the appearance of asymp- 
totic unitarity, as well as the study of the distance 
\\p{t) — p{t)\\, has a major advantage with respect to 
those based on p (t). Namely it allows to discriminate be- 
tween approximate unitary evolution and the usual reach 
of an asymptotic state for which p = 0. 

We would like to end this section by stressing 
the (almost obvious) relation of asymptotic unitarity 
with the quantum-information concept of noiseless or 
decoherence-free subspace/systenP^. The quantum net- 
works considered in this paper are of hybrid type, namely 
some of the inter-site couplings are coherent i.e., hop- 
ping, and other are incoherent i.e., irreversible transfer 
described by L. On the other hand, the dynamics re- 
stricted to the range of the projection P :— li2^ is 
unitary because, as noticed in the above, Cl{PpP) = 0. 
This means that the range of P is indeed a decoherence- 
free subspace. Now the dynamics is such that, for ap- 
propriate initial conditions limt_i.oo 71,3 (t) = 1 or equiva- 
lently limt_j.oo \\Pp{t)P — p{t)\\ — . This means that the 
asymptotic state belongs to the range of P, which in turn 
implies the unitary nature of the long-time dynamics. 



V. APPLICATIONS TO LHl-RC COMPLEXES 

In this section we want to check if and how the ef- 
fects studied so far can survive in more realistic networks. 
Specifically, we will consider models which can be rele- 
vant for the description of energy transfer in photosyn- 
thetic systems. 
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Figure 7: A'' particles on the ring interact via dimerized near- 
est neighbor hopping constants — t (^1 + S (—1)'^ . Each 

of these particle can hop to the central site (A^ + 1th) with 
hopping constant J. The central site transfers excitations in- 
coherently to the reaction center via L — -v/70"jv+i''"flc- "^^^ 
reaction center itself is connected to a (2s +1)— dimensional 
"battery" via Lb — x/Tt^^-RC^baf '^^^ effect of external de- 
grees of freedom is schematized by incoherent dissipation and 
dephasing processes (thin, blue, wavy lines). In actual LHl 
complexes the sites on the ring are bacteriochlorophylls, and 
N = 32. 



A. Congestion effect 

Our motivation for the study of the "congestion effect" 
originated from a careful analysis of the structure of the 
reaction center in LHl-RC complexes. In most photosyn- 
thetic bacteria, photons are captured by light-harvesting 
antennae where a particle-hole exciton is created and car- 
ried to the reaction center (RC) where eventually a redox 
reaction takes plac^. In the LHl-RC complexes present 
in purple bacterigPS the light harvesting complex and the 
RC form a compact core unit. Typical transfer times of 
excitations to the RC are of the order of hundreds of pi- 
coseconds. A cartoon picture of the LHl-RC complex 
is shown in figure [7] Yellow spheres represents the bac- 
teriochlorophylls forming LHl. In the purple bacterium 
Rodobacter sphaeroides, there are 32 bacteriochlorophylls 
(BChl) displaced on a ring surrounding the reaction cen- 
ter. In figure [7] we display a possible structure for the 
RC. Instead of treating the RC as a simple two- level sys- 
tem, as typically done in the literature, we replace the 
RC with a structure containing two qubits and a d-level 
system which we call a "battery". In purple bacteria 
this structure has to be imagined sitting at the center of 
the ring. The first of these qubits (the N + 1th) inter- 
acts via coherent dipole-dipole hopping with the BChls 
of the ring. Excitations are then transferred at a rate 7 
to what we call reaction center. In turn, the RC itself 
is connected to larger c?-level system (d = 3 in our sim- 
ulations) via irreversible transfer at a rate jt- It is the 
interplay between the two timescales 7"^ and and 
their relation to the transfer efficiency, that we want to 
analyze here. 

The master equation for the whole system is of Lind- 



blad type: p = —i [H, p] + £tot (p)- For what we said so 
far, the incoherent part is given by £tot = Cl + J~-Li, + 

^noiso with L = y/j'^N+l'^RC ^^^^ Lfc = Vlb<^RC^baf 

Dissipation and dephasing effects are taken into account 
via incoherent terms acting on the sites of the ring 



N 



with ij,diss = VTdil 



-1 ^j.diss ' ^3 

and ij,dcph = -y/7dcphi^j • 

Regarding the Hamiltonian of the ring degrees of free- 
dom, w e refe rred to the detailed structure of couplings 
given i] j35|36 | ^Y^e most salient feature emerging from 
the data oP^ is that the couplings present a dimerized 
structure: strong coupling t+ = t{\ + 5) alternate with 
weak ones t- =t(l — 5). Indeed, instead of using all the 
couplings tij reportecP^, almost the same band struc- 
ture can be obtained using only a nearest neighbor de- 
scription with a dimerization of 5 = 0.12. Some groups 
have suggested the possibility that dimerization might 
favor the transfer efficiency'^''. Our choice of resorting to 
a dimerized nearest neighbor hopping structure has the 
additional advantage of making the system scalable to 
different sizes TV. Hence our choice for the Hamiltonian 



N 

This represents N particles on a ring hopping between 

neighboring sites with constants tj = t (l + 5{-iy^ and 

to a central site + 1 with hopping constant J. We will 
also add static random diagonal noise {H — )■ H+'Y^- e^rij) 
to inhibit the possible appearance of decoherence-free 
subspaces which can limit the efficiency of transfeJ^. 

The results of our simulations are shown in Fig. [Sj We 
initialize the system by starting with a pure Dicke state 

for the ring while keeping all other sites empty. This 

, ^ , — 1/2 

means the initial state is l-fAo) = ( I (o'tot)" |0) 



where cr^^j — X^jLi '^^ refers only to the ring sites and 
|0) is the empty state for the whole system. The choice 
of an initial Dicke state is natural for a series of reasons. 
First it allows to treat initial states with general defi- 
nite particle number n < N . Second, Dicke states are 
symmetric under permutation, thus carrying no net mo- 
mentum. If the photon's wavelength is larger than the 
size of the LHl complex, the excitations created must be 
a completely delocalized k = Q object. In any case, since 
only the k = Q component of the ring couples to the cen- 
tral N -t- 1th site, transfer in the antisymmetric channel 
k = TT, being a higher order process, is highly suppressed 
and gives much lower transfer efficiency^. 

We first performed simulations on a "clean" system, 
i.e. with no dissipation or dephasing present. In Fig. [8] we 
plotted the population of the reaction center (normally 
called efficiency 77 in the literature) as a function of time 
for different values of "fb- Looking at the upper panels 
of Fig. [8j the situation is completely analogous to the 
congestion effect observed in our simple toy model (see 
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Figure 8: "congestion effect" in light-harvesting complexes. 
Upper panels: clean system, no dissipation and dephasing. 
The ring has N = A sites, hopping constants are t — J = 
l{meV) and dimerization is 5 = 0.12. Diagonal static noise 
of the form e„ = t cos (en) is added. The RC transfer rate 
is set to 7 = 0.3ps~^. Lower panels: same parameters plus 
dissipation and dephasing 7diss = 7deph ~ 0.03ps~^. Left 
panels: the initial state is a two-particle Dicke state for the 
ring, other sites are empty. Right panels: the initial state is a 
three-particle Dicke state for the ring, other sites are empty. 
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Figure 9: A'^ + 2 = 6 sites. On the upper panels we plot the 
populations of the ring and of the RC as a function of time 
(arbitrary units). Bottom panels: parametric plot of the real 
and imaginary parts of a matrix element {tl),n\{l\p (t) 
for certain m, n. | -(/'„) are the Hamiltonian eigenstates. Left 
panels: the Hamiltonian has t — J — 1, S — 0.12, diagonal 
static noise ep — t cos (ep) , and no dissipation or dephasing. 
Excitations are transferred to the RC at a rate 7 = 0.2. Right 
panels: same parameters, but on the particles of the ring acts 
dissipation and dephasing with 7deph = 7diss = 0.01. 



Figures |2j [3]). As long as we start with a number of 
excitations which can be accommodated in the battery, 
they will all flow to the battery for 7fc ^ (left panel). 
When we start with 3 particles in the ring we see again 
the appearance of a non-monotonic behavior between 
and 7 which shows up as a valley at large times and 
lb 1^ 1 {lb smaller than, but of the order of 7). When we 
add additional decoherence in the form of dissipation and 
dephasing the situation is only quantitatively changed. 
The valley due to the "congestion effect", although less 
pronounced, is still visible in the bottom right panel of 

Fig. ID 



B. Asymptotic unitarity 

To study asymptotic unitarity the "battery" is an un- 
necessary complication. Therefore, we will use the same 
model of Fig. [7] without the battery site and the corre- 
sponding jump operator. This leads to a network of A'^-l-2 
qubits where the N + 1th site is connected to the RC via 
irreversible transfer at a rate 7. As done previously, we 
will use an n-particle Dicke state as the initial state. Let 
us first consider the case where the only incoherent term 
is the one transferring particles from the central site to 
the RC. In this case the dynamics becomes exactly uni- 
tary when the RC is full. Simulations on a network with 
N + 2 — Q qubits are shown in Fig. |9] We also show the 
effect of dissipation and dephasing, though one order of 
magnitude smaller than the RC transfer. For short times 
the evolution is the same as for the clean (i.e. no dissipa- 



tion and dephasing) case, however for time of order 7^^^, 
dissipation sets in and the parametric plot for a generic 
matrix element {n\p{t) |to), spirals down to zero (Fig. [9] 
bottom right plot). 

The conclusion of this section is as simple as it is in- 
triguing, in view of potential applications to biological 
systems. If the time-scale 7jjgg is large enough, there may 
exist a time window T^oiax < t < l^^l^ in which quantum 
effects are not only visible but the dynamics is effectively 
unitary! In our models T^d^x is the time needed for the 
RC to get filled, and is of the order of Tj-oiax ~ 7~^. Even 
more important is the fact that the incoherent transfer 
to the RC must shut down when the RC is full. Consid- 
ering the LHl-RC complex, the sep aratio n of time-scales 
does indeed occur. For instance irpSEIl the dissipation 
is four orders of magnitude smaller than the RC charge- 
separation rate. Whether the RC shuts down when it is 
occupied, although plausible, is much harder to assess. 



C. Staircase effect 

Here we want to show that the staircase effect, studied 
in Sec. IIIIl survives in more elaborate networks. We will 



study this effect in the model depicted in Fig. 10 The 
model consists of an open chain of N sites hopping co- 
herently between nearest neighbors, i.e. the Hamiltonian 
is 



N-l 
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in 



= ^7dcphnj as defined 



1- ^- ^- ) 
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Figure 10: sites interact via a nearest-neighbor hopping 
Hamiltonian. Particles are injected, and respectively expelled 
incoherently at rates 7in, 7out on the first and last sites . On 
top of this basic structure we can add static diagonal disorder 
and dissipation as well as dephasing (symbohzed by blue wavy 



arrows) L-a 
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Figure 11: J = 1 and 7in = 0.2, 7out =0.3. First row panels: 
neither dissipation nor dephasing, and no static noise. Second 
row panels: addition of diagonal static noise e„ = J cos (en). 
Third row panels: static noise plus dissipation and dephasing 
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Particles are injected into the first site of the chain via a 
jump operator Ljn = ^JyirS^X and taken away at the last 
site via Lout — VTout . On top of this basic framework 
we add different layers of complexity. First we can add 
some static random diagonal noise, i.e. we add site de- 
pendent energies to the coherent part H ^ H + J2j ^j^j- 
Second we can also include dissipation and dephasing 
acting on the inner sites of the chain by adding the fol- 
lowing superoperator: £noiso = EjL^^ ^Lj^^i.s + ^L,_d.ph 



{Lj,diss = y/ldi^cr~ and L 
previously) . 

The picture that we have is the following. Through the 
coherent part of the evolution, excitations travel in the 
chain in packets of quasiparticles at velocities Vk = O ( J) 
(/c is a quasi-momentum label). This introduces a lag 
timescale Tq ^ L/v ^ O {L/ J), which is the time needed 
for an excitation to travel from one side of the chain to 
the other. From Fig.[ll](all panels) we see that, when the 
population at the injection site increases, the population 
at the expulsion site stays constant during this time-lag 
To and vice versa. Considering the central and lower 
panel of Fig. [TT] we can appreciate how robust the effect 
is with respect to various type of "perturbations". The 
addition of static random noise has the effect of localiz- 
ing states and shuffling the single-particle dispersion e^. 
Both of these effects destroy the picture of wavepackets 
traveling at constant velocity, in that both the traveling 
times and the dispersion of the wave-packets increase. 
Instead the addition of dissipation (and dephasing) to 
the network, mostly has the effect of relaxing the system 
at a faster rate. As long as the system has not relaxed 
the effect remains visible. Comparatively, the presence 
of static coherent noise hinders the stair-case effect more 
to dissipation and dephasing. 



VI. CONCLUSIONS 

Inspired by the models which are recently being used to 
describe energy transfer in photosynthetic pigments, we 
have identified and discussed a few effects arising in quan- 
tum networks with coherent (Hamiltonian) as well as in- 
coherent (Lindblad) coupling between the nodes. For the 
reader's sake we summarize here below these basic effects 

1. Congestion effect. The incoherent transfer of exci- 
tations is inversely proportional to the population 
in the reaction center. This is due to the hard-core 
nature of the excitations that effectively reduces 
the amplitude of the jump operator as the reaction 
center fills. 

2. Asymptotic unitarity. Coherent, unitary evolu- 
tion may emerge out of a dissipative, incoherent 
dynamics. This happens if states which annihilate 
the incoherent part of the dynamics can be reached 
during the time evolution. For this effect to be 
observable one needs a separation of time-scales, 
Trciax *C Tdiss • Such Separation of time-scales does 
take place in some photosynthetic systems e.g. in 
the LHl-RC complexes present in purple bacteria. 

3. Staircase effect. This effect refers to a situation 
in which particles are injected incoherently, travel 
coherently along a given chain, and then are ex- 
pelled (or digested) at a certain rate at the other 
end of the chain. The effect of the coherent part is 
to introduce a time-scale Tq — O (L/v) = O (L/J), 
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(L is the system size, v the velocity of excitations 
and J is the energy scale of the coherent network) . 
Tq is roughly the time needed for the excitations 
to travel from one side of the chain to the other. 
The peculiar feature emerging from the dynamic 
evolution, is that when the population at the injec- 
tion site increases, the population at the expulsion 
site stays constant during this time-lag Tq and vice 
versa. This effect results in a step-like behavior 
in the parametric plot of the injection/extraction 
populations. 

The effects we analyzed in this paper can be traced back 

to very simple mechanisms displayed even by networks 
composed by only few qubits. We provided analytical 



solutions for these toy models and showed numerical ev- 
idence that these effects survive in more elaborated net- 
work such as those modeling energy transfer in purple 
bacteria. Clearly, further investigations are in order to 
establish the relevance of the elementary calculations pre- 
sented in this paper to the newborn field of quantum 
biology. 
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